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Abstract: We consider the Gibbs-measures of continuous-valued height configurations on the 
(i-dimensional integer lattice in the presence a weakly disordered potential. The potential is 
composed of Gaussians having random location and random depth; it becomes periodic under 
shift of the interface perpendicular to the base-plane for zero disorder. We prove that there 
exist localized interfaces with probability one in dimensions d > 3 + 1, in a 'low-temperature' 
regime. The proof extends the method of continuous-to-discrete single-site coarse graining that 
was previously applied by the author for a double-well potential to the case of a non-compact 
image space. This allows to utilize parts of the renormalization group analysis developed for the 
treatment of a contour representation of a related integer- valued SOS-model in [BoKl]. We show 
that, for a.e. fixed realization of the disorder, the infinite volume Gibbs measures then have a 
representation as superpositions of massive Gaussian fields with centerings that are distributed 
according to the infinite volume Gibbs measures of the disordered integer-valued SOS-model 
with exponentially decaying interactions. 
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I. Introduction 

The study of interface models from statistical mechanics, continuous as well as discrete 
ones, with respect to their localization vs. fluctuation properties, is an interesting topic in 
probability theory. In this paper we study the problem of continuous SOS-interfaces in random 
potentials that are random perturbations of periodic ones and prove stability of the interface 
in dimensions d > 3 + 1 (as suggested by the heuristic Imry-Ma argument, known for long to 
theoretical physicists). 

A related stability result has been proved before for the simpler discrete version of such a 
model with nearest neighbor interactions in [BoKl]. The proof uses a renormalization group (or 
spatial coarse-graining) procedure that was based on the technique of Bricmont and Kupiainen 
that was developed for the Random Field Ising Model [BK]. The issue of this note is thus to 
clarify what to do with additional (possibly destabilizing) fluctuations of the continuous degrees 
of freedom. 

An analogous problem was investigated in a recent paper by the author [K4] in the simpler 
case of a random double-well potential, where ferromagnetic ordering was shown in d > 3 
(under suitable 'low temperature' and 'weak anharmonicity' assumptions on the potential). 
The key point here is to construct a suitable stochastic mapping from continuous to discrete 
configurations and study the image measures under this mapping. In the double-well case this 
mapping is just a smoothed sign-field indicating what minimum the continuous spin is close 
to. The image measure could then be shown to be an Ising-measure for a suitable absolutely 
summable Hamiltonian. It can be controlled by the known renormalization group method of 
[BK] . This is clearly in favor of running a suitably devised renormalization group transformation 
on the (in this context unpleasantly rich) space of continuous configurations. 

The purpose of the present paper is to study the difficulties of infinitely many minima in the 
potential. The stochastic mapping we will apply to the continuous spins will now be a mapping 
to integer- valued spin configurations. As opposed to [K4] the mapping will now also depend on 
the realization of the disorder. 

To explain the method in the simplest non-trivial context, we have decided to choose a 
specific potential that is the log of sums of Gaussians. The treatment of this potential provides 
the basic building block of the analysis also for more general potentials in that it explains the 
occurrence of the phase transition and the structure of the contour models that will arise. It 
corresponds to having vanishing 'anharmonic corrections'; how those anharmonicities (that are 
present for more general potentials) can be treated by additional expansions is explained in 
detail for the double- well case in [K4], so that combining those methods with the ones from the 
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present paper should yield stability for a larger class of continuous interface models. 

This restriction also allows us to obtain particularly nice 'factorization-formulas' for the 
continuous-spin Gibbs-measures in finite and infinite volume. They have some probabilistic 
appeal and clarify the structure of the coarse graining transformation we use. In particular 
we can describe the infinite volume Gibbs measures in terms of the 'explicit' building blocks 
of random discrete height measures and well-understood (random) massive Gaussian fields (see 
[1.7]). 

Here is the model. We are aiming to investigate the Gibbs measures on the state space 
IR Z of the continuous spin model given by the Hamiltonians in finite volumes KdZ d 
E ~ A(mA) = | ^ {m x -m v f + \ £ (m x -m y ) 2 + Y,V x (rn x ) 

{i,s)CA a:GA;y€dA x£A 

d(x,y) = l d(x,y) = l 

for a configuration € IR A with boundary condition nig a- We write <9A = {x G A c ; 3y G A : 
d(x,y) = 1} for the outer boundary of a set A where d(x,y) = \\x — y\\i is the 1-norm on IR d . 
q>0 will be small. 

The random potential we consider is given by 
V x {m x ) = - log 
m* x (l)=m* {l + d x (l)) 

The disorder is modelled by the random variables (Vx(h)) x€Z d. h€Z and (d x (h)) x€Z d. h€Z , de- 
scribing the random depths of the Gaussians and the random deviations of the centerings of 
the Gaussians from the lattice m*/Z. The unperturbed potential thus takes its minima for 
m £ m*/Z, the fixed parameter m* > being its period. Later it will have to be large enough. 
(Note that the curvature of the potential is of the order unity for large m* ; thus the curvature 
really has to be large on the rescaled lattice where the potential has period 1.) 

We will simply take the {Vx{h)) x€Z d. h€Z as i.i.d. random variables with distribution IP V 
and the (d x (h)) x€Z d. heZ as i.i.d. random variables with distribution IP^, independent of the 
r/'s. Furthermore we impose the smallness conditions 

_ * 2 

(i) IP[\vAh)\ >t]<e (ii) \ Vx (h)\ < 5 V 
tj_ 

(iii) IP [d x {hf >t]<e 2 °d, ( iv ) \d x (h)\ <5 d <\ 

where cr^,(j^ > will be sufficiently small. 

An assumption of the type (iv) is natural, since it just states that the shifted wells stay 
away from each other and don't merge or even cross. The assumption (iii) is less natural (and 
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l(m x -m*(l)) 2 +ri x (l) 



where 



(1.2) 
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not really essential). Moreover we will need in the proof that Sd, S v be sufficiently small, which 
is just to simplify the structure of the contour representation we will derive later and could be 
bypassed, see below. 

To make explicit the local dependence of various quantities on the disorder variables we 
write lo x = (d x (h),r) x (h)) heZ for the 'disorder variables at site x' and put uj a = {oo x ) xeA - 

A more general setting could of course be to consider V x that are stationary w.r.t. a discrete 
shift in the height-direction and satisfy some mixing condition. Also the i.i.d. assumption in x 
could be weakened. 

We use the following notation for the objects of interest, the finite volume Gibbs-measures 
/j,™ 9A,1jJA , defined in terms of their expectations: 

^ A (f) = l^r A [ dm A f(m A ,m A c)e- E ^ A ^ where 

JmA (1.3) 
A = / dm A e- E ^ A ^ 

JlR A 

for any bounded continuous / on IR Z (continuity is meant w.r.t. product topology). Most of 
the times we will put zero boundary conditions m x = (for all x € % d ), writing simply l^ A U3A - 
(Due to stationarity that's the same as putting m x = m*l for any fixed I € Z.) 

Then we have the stability result 

Theorem 1: Let d > 3 and assume the conditions (i)-(iv) on the disorder variables. Then, 
there exist qo > (small enough), So > (small enough), To < oo (large enough), Oq > (small 
enough) such that, whenever S V ,S^ < So, q(m*) 2 > To, and of ff := u d + a 2 < ctq, the following 
is true. 



There exists an infinite volume random Gibbs-measure ^ that can be obtained as the weak 
limit /x w = limjvfoo ^\^ N along a non-random sequence of cubes A^r. The measure describes a 
continuous interface localized around the base plane; its 'roughness' is bounded by 

TE^ (m 2 Xo ) < 1 + (m*) 2 ^ e - const ^ + e'^^j (1.4) 
for any x , with (3 = const x min < log ^,qm* 2 ^ lo ^» ^ \ and an exponent k > 0. 



Remark: So, measured on the scale of the period m*, the roughness is in fact a very small 
number. The term 1 has to be present since it describes the true fluctuations of a continuous spin 
in an individual well of V x . The quantity q{m*) 2 gives the true order of magnitude of the minimal 
energetic contribution of a pair of nearest neighbor heights in neighboring potential wells, so it 
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can be viewed as some basic temperature variable that has to be large enough. It appears in the 
definition of $, with some minor logarithmic deterioration that we need for technical reasons. 
The log ^-contribution comes from a high-temperature expansion, to be explained later. 

Essential for the analysis is the following local transition kernel T x ( ■ I ■ ) describing a single 
site coarse graining from a continuous height m x 6 IR to an integer height h x £ Z. It is defined 
by keeping the Gaussian for I = h x , 

T%- (h x \m x ) = — , (1.5) 

1 ' Worm. (m x ) 

where Norm.(m x ) is chosen to make T x x (h x \m x ) a probability measure on /Z, for fixed m x . Note 
that this object depends on the disorder through oj x (as opposed to our proceeding in the double- 
well case where a simpler analogous kernel was non-random). So, for fixed m x , the random 
probability weights (T x {h\ m x)) h€Z ar e integer samples of a randomly perturbed Gaussian. The 
reader may also note the complementary fact that, for fixed h x , the probability T x (h x \m x ) is 
bounded by Gaussians from below and above and has a unique absolute maximum as a function of 
m x . The maximizer is close to m*h x . [See Lemma 2.1 for these statements]. Now, the simplicity 
of our basic log-sum-of-Gaussian potential (1.2) lies in the fact that Norm.(m x ) = e~ Vsc ^ m ^\ [For 
more general potentials V x (m x ) this equality will acquire error terms ('anharmonicities') that 
lead to additional expansions.] We use the same symbol, T ^dh\\mjC) = EI^eA T x x (h x \m x ), for 
the kernel from IR A to /Z A , and also in the infinite volume. 

Before we put down more results in a precise way, let us describe in an informal way in 
symbolic notation what we are about to do. Starting from a (finite volume) Gibbs measure jj,{dm) 
we look at the joint distribution M(dh,dm) := T{dh\m)ix(dm) on integer heights configurations 
h and continuous heights m. Then we analyse \i with the use of Bayes' formula: We have 
/j,(dm) := J u(dh)M{dm\h) where u(dh) = J n{dm)T(dh\m) is the ^-marginal. It is of course 
a completely general (and a-priori empty) idea to look at distributions in suitably extended 
space that can only be useful for natural choices of this space. In our case we succeed with the 
control of v{dh) since we can obtain a contour representation that can be treated by the spatial 
renormalization group. The conditional probability M(dm\h) is nice for the specific choice of 
the potential; it is just a Gaussian distribution. (M(dm\h) would be more complicated for 
perturbations of the potential, but the above decomposition would still be a successful one.) 
Our results, to be described below, will then concern the approach of the thermodynamic limit 
of v. We will also have to clarify the interplay of the thermodynamic limit with the above 
formulas. 

The following theorem describes how the control of Gibbs-measures on the integer heights 
carries over to the control of the Gibbs-measures on the continuous heights, under some harmless 
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additional condition. 



Theorem 2: Suppose that the discrete height measures z/£ := T \jJ*™ aA ' WA j converge locally 
along a sequence of cubes An, centered at the origin, to a limiting measure u^{dh Z d) for a 
sequence of boundary conditions which is uniformly bounded, i.e. we have snp xeZ d. N \fhQ\ N - x \ < 
M , for some M < oo. (That is, convergence takes place for expectations of all bounded local 
observables.) Assume moreover that we have the site-wise summability 



S ^ P Yl i 1 - 1 A Z d ) x ^A N \ h y \ =■ K xH < OO 



2/GAj, 



(1.6) 



for a sequence of increasing cubes An, for all x € /Z d for a.e. configurations of the disorder u. 
Then the measures fj,™ 9A,u ' A converge locally to the infinite linear combination of Gauss measures 
given by 



M w := Jif(dh z *)tf[{l-qA x *) l m* zd {h zd )-(l-qA zd ) 1 



(1.7) 



That is, convergence takes place for sequences of expectations of bounded measurable f(mv) that 
depend only on spins in the finite volume V. 

The symbol M a; (1 — </A z .i )~ ' denotes the massive Gaussian field on the infinite lattice 
Z d , centered at a G lR zd with covariance matrix given by the second argument (so that we 
have e.g. JM a; (1 - gA Z d) _1 (dm)(m x - a x )(m y - a y ) = (1-gA^d)" 1 ). 

J L J x,y 

Remark: Note that the random quantities K x (lo) will typically not be bounded uniformly 
in x. In fact, even a localized interface will have unbounded fluctuations around regions of 
exceptionally large fluctuations when considered in the infinite lattice. Of course, (1.6) is implied 
by swp Nty IEi>%\hy \ < oo. 

Remark: We stress that Theorem II does not only apply to the 'flat' interfaces that we 
investigate here but also to more 'exotic' Gibbs-measures. So, e.g. the (supposed) existence 
of Dobrushin-type integer-height Gibbs-measures (that are perturbations of a flat interface at 
height in one half-space and a flat interface at height H in the complement) would imply the 
existence of corresponding continuous spin Gibbs-measures. 



The organization of the paper is as follows. In Chapter II we prove the 'factorization' 
Theorem 2, starting from its finite volume version Lemma 2.1. In Chapter III we derive the 
contour representation of the integer height model (see Proposition 1), starting from the finite 
volume Hamiltonian (3.1). In Chapter IV we conclude to prove Theorem I from these results 
applying the spatial renormalization group construction from [BoKl], [Kl] on the contour model 
representation. 
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l+2<5 d 
l-2S d 



+ 25 v > . We have that 



II. The Joint distributions of continuous and integer heights 

Before we get started, let us make explicit some (concentration-) properties of the random 
transition kernel to get some intuition for it. The elementary proof is given at the end of the 
chapter. 

Lemma 2.1: For any realization of the disorder satisfying the bounds (ii) and (iv) below (1.2) 
the fixed-h x (random) probability m x — > T x x (h x \m x ) = 2 ^ \ \ t) has a unique 

absolute maximum. It has no other local maxima. The maximizer lies in the symmetric interval 
about m*h x with radius ^jjT^y^ + m *2{i-S d ) j lo g 

const e -5(™*-<(M) 2 < Tx (h x \m x ) < Const e -5( m --<(M) 2 (2.1) 
with const > 0, Const > depending only on a, m* , 5d, S v . 

Now, the simplicity of our choice of the log-sum-of-Gaussian potential lies in the fact that 
the joint distribution on continuous heights and integer height can be written in the form 

M ~(dm A ) [J T x {h x \m x ) = -J^e- H ^ A ' hA ^dra A (2 . 2) 

xeA Z A 

where 

= ^ E i m x- m y) 2 + 7} ( m x-rh y f + ^J2(m x -m* x (h x )) 2 -J2n x (h x ) ( 2 - 3 ) 

{x,y}CV x€V;yGOV igV X^V 

d(x,y) = l d(x,y) = l 

is quadratic in m, for fixed h. This is due to the cancellation of the normalization in the 
transition kernel against the exponential of the potential. We remark that, for potentials that 
can be viewed as perturbations of our specific log-sum-of-Gaussians the formula would acquire 
error terms and the present formula is the main contribution of a further expansion. 

We will now rewrite the joint distribution as a product of the marginal in the integer 
heights and the conditional distribution of the continuous heights given the h. We see that the 
m-distribution conditioned on a fixed value of h is Gaussian. The /i-marginals on the other hand 
can be computed by a Gaussian integration over m\: Since the quadratic terms of the above 
integral are /i-independent this Gaussian integration yields 

I u-»"SA'"A.' l A / -> ;„( . rr"''8A> u 'A>' l A/"„ 1 \ 

/ rfmAe^A 3 * A W =CAXe -< a „Aif A (m A ) ( 24 ) 

Jm A 

with a constant C\ that does not depend on h\ (and uj\). 
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By multiplying and dividing the r.h.s. of (2.1) by (2.3) we get after a little rewriting of the 
Gaussian density, conditional on the h: 

Lemma 2.2: The finite volume joint distribution of continuous and integer heights can be 
written as 

^>" A (dm A ) J] T x (h x \m x ) = u^(h A ) M \mt aMhA ; (l - gAf)" 1 ] (dm A ) where 



u^(h A )=T^^) (h A ) 



e 



inf mA€IRA H^-- h \ mA ) (2.5) 



with the random centering 

m^ A ' hA = (1 - ^Aa)" 1 (m* A (h A ) + qd A , dA m dA ) (2.6) 



Here d At Q A is the matrix having entries (d Aj d A ) XjV = 1 iff x G A and y G <9A are nearest 
neighbors and zero otherwise. Aa is the Lattice Laplacian (with Dirichlet boundary condition) 
in the volume A. 

Side-remark: We like to point out that the enlargement of the probability space by the 
introduction of auxiliary integration variables and conditioning on the latter ones can be found 
in various places in statistical mechanics: 

1) In the renormalization group analysis one studies the Gibbs-distributions fiiT) of the 
variables V of the system with the aid of a mapping to spatially coarse-grained variables V 
by means of a transition kernel T{dT'\T). This gives rise to a joint distribution M(dT,dT') = 
n(dT)T(dT'\T). Reversing the order of conditioning gives M(dT, dY') = ^(r")M(dr|r'), the idea 
being that the 'renormalized' measures v(T') are easier to study than the measures /u(r). (Of 
coarse, even if this is true, the conditional distribution M(dT\T') must also be controlled.) 

2) The introduction of artificial integration variables is a commonly used trick also for the 
analysis in quadratic mean-field models (known here as Hubbard-Stratonovitch transformation) . 
In fact, the analogue of formulas (2.5) looks as follows for the simplest candidate, the usual mean- 
field Ising ferromagnet. Its Gibbs distribution on the spins (<Ti)i=i,...,N £ {—1, 1}^ =: is given 
by Hn(o-) = e^ mjv(CT ) 2 /Norm. where m N (a) = YliLi a i- Consider the variables (<r, rh) in the 
enlarged probability space VLxIR (with a new smoothed out magnetization variable rh) that are 
distributed according to the joint distribution M(a,dm) = e f 3Nm (< 7 )- m e -^ m ' 2 dm/Norm.. Then 
the marginal distribution on the a is the desired Gibbs-distribution; in fact we have M(a, dm) = 
n(a)T(m\a) with T(dm\a) = e~^~^ m ~ m ^ 2 dm/Norm, describing the smoothed out averaging 



8 



over the whole lattice. To analyse the Gibbs-distribution the joint distribution is then written 
by reversing the order of conditioning in the form M(a,drh) = v(dm)M{a\rh) where v(drh) is 
the marginal on the m and M(a\m) = M^a^rh) where M^a^rh) = e /3mo-;-iog2cosh(/3m)_ The 
latter kernel is clearly trivial, the distribution u(drh) is treated by a saddle-point method. 

Likewise, our strategy in the present problem will now be to control the ^-distribution 
in the thermodynamic limit and get the Gibbs-measures of the continuous spins by summing 
(2.5) over h. Assuming this control over the integer heights we must however also control what 
happens to (the /i-average over) equation (2.5) under the thermodynamic limit if we apply it to a 
local function /(my), depending on continuous heights m x only for x G V, V being a fixed finite 
volume. Note that the Gaussian describing the conditional distribution of the continuous heights, 
given the integer heights, has some A- dependence both through its centering and the covariance 
matrix. Further, its dependence on h is not finite range (albeit strongly decaying unless the 
integer- heights are getting very large.) So we need some extra condition on the convergence 
of z^A-measure and a little work to deduce the implication of the desired convergence of the 
/iA-measure. 

The precise result of this is given in Theorem 2 (see Introduction), that we are going to 
prove now. While doing so, we will also prove the following 

Addition to Theorem 2: Assume the site-wise existence of all exponential moments 



=: K x {cj,s) < oo, (2.7) 



Then the convergence limAfoo ti™ 9A '^ A (/v) = h^ifv) takes place also for all local observables 
/(my) that do not increase faster than exponentially; i.e. there exists a constant A > s.t. 
\f(m v )\ < e A ll mv " 2 . 

Proof of Theorem 2: We start with the proof of the theorem under the 'site-wise 
summability assumption' (1.5). We must control large realizations of the h's (that are however 
improbable w.r.t is, under this assumption.) Let / denote any measurable function of my, we 
assume for simplicity that / is uniformly bounded by 1. To produce a local observable (of the 
integer heights) we cut off the long range dependence of the Gaussians on the integer height h 



sup v% 

A 



d) \h y \ 
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outside some volume A 2 that satisfies VcA 2 cA. We use an e/3-trick to decompose 



m 



(1 - qA z .i) 1 m* zd (h Z d) ; (1 - qA Z d) 1 (/) 



AT 



(1 - gA^)" 1 (m* A2 (h A2 ) , A c) ; (1 - qA^y 1 } (/) 



+ 



K - [(1 - qA z ,r y (m* A2 (h A2 ) , A c) ; (1 - qA^y 1 } (/) 



( 1 - <lA z ,i ) 1 m* d (h Z d);(l - qA Z d) 1 (/) 



-A/" 



(1 - «yA z ,)-' (m A2 (h A2 ) , A c) ; (1 - gA^)" 1 ] (/) 



(2.8) 



We will show that the r.h.s. can be made arbitrarily small be choosing at first the auxiliary A 2 
and then A large enough. Indeed, the middle term on the r.h.s. converges to zero with A j oo, 
for any fixed A 2 , due to the assumption of weak convergence of the The remaining task is 
to control the two error-terms; for this we need the condition (1.5) [resp. (2.7)]. 

We look at the first term on the r.h.s. more carefully. The last term is treated in a 
similar fashion. We need some continuity properties of |V (-dimensional Gaussian expectations 
considered as functions of their means and covariances. The following estimate will do, both for 
bounded observables and observables that are only exponentially bounded. 

Lemma 2.3: Let J\f[a, E], jV[a',S'] denote two \V {-dimensional non- degenerate Gaussians 
with mean a, a' € IR V and covariance £, £' € IR VxV . Assume that /(my) is an observable that 
doesn't increase faster than exponentially, i.e. \f{m v )\ < e A H mv l' 2 for some A > 0. Then we 
have the following estimate 



M [a; £] (dm v )f(m v ) - J M [a; £'] (dm v )f(m v ) 
<2 |vVNI 2e ^ 



det£ 
det £' 



+ 



det£ 
det £' 



((25 + ||a|| 2 + ||a'|| 2 ) HE" 1 ^ ||a - a'\\ 2 + 2 (S 2 + Ha'Hl) x ||E" 



+ 2 |y| e AS- < 5 -^' +2 \v\ 



AS- 



(S-lla'lb) 2 



(2.9) 



where g(x) = xe x , for any S > max{||a|| 2 + ATrE, ||a'|| 2 + ATrE'}. 

Proof: To show the Lemma we decompose the range of integration into a ball of radius S 
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and its complement; to control the corresponding integrals we need a simple application of the 
exponential Markov inequality in the form of 



Lemma 2.4: Let J\f[a, £] denote a Gaussian with mean a £ IR V and covariance £ £ 
(i.e. J M[a, T,](dmv) (m x — a x ) (m y — a y ) = T, x ^ y ). Then we have 

N [a; £] (dm y )e A " m ^ 2 < ^^e^e^ and 

.\'[a;E](dm y )e A|!mv||2 l|| mv ,|| 2 > 5 < ^e^'^^ 1 for S > \\a\\ 2 + ATr£ 
For S = II a II 2 + ATrE the two r.h.s. coincide. 



J- 



(2.10) 



Proof: Write J'Af[a;T,] (dm v )e x ^ mv ^ 2 < e A H a " 2 / [a; E] (dmyje*""^ -0 " 3 and denote by of 
the eigenvalues of E and by e, the corresponding eigenvectors. Then we may write 



f M [a; E] {dm v )e x W mv - a ^ < [ [J I M [0; of] (dm, 

J \i=l,...,\V\ J 



NgAlm.l < 2 |V| e 2 



V „ 2 



(2.11) 



which proves the first estimate. The second one is a corollary: Write, for nonnegative Ai, 



J M [a; E] (d 



my)e A|| m v|| 21 < g - 



(A + A^TrS 



Kl5 y TV [a; E] (dm y ; 



,(A+Ai)||m v ||a- 



; ' my ||,>.S 



(2.12) 



< e -AiS 2 |V| e (A+A 1 )||a|| 2e 

where the last inequality follows from the first claim. Minimizing the r.h.s. yields the claim (the 
optimal value of Ai being Ai = lS '~[^^ 2 — A). The requirement that Ai be positive leads to the 
given range of allowed S. C" 



We continue with the proof of the Lemma 2.3 writing 

J M [a; E] (dm v )f(m v ) - J N [a'; £'] (dm v )f(m v ) 

J Af[a;X}- J Af[a';X'] \ (dm v )f(m v )l ]mv] < s 



(2.13) 



+ 



j A/" [a; E] . 



,a|!«h 



'1 



\my \>S 



+ |;V[a';S']e A ll m -lll| mv |> 5 



The last two terms are estimated with the help of the above lemma, leading to the last line of 
(2.9). The first term can be estimated simply in terms of differences of the Gaussian densities: 



dm\ 



< sup 

\m v \<S 



e 2 



i<(mv— a),S 1 (mv— a)> g — \ <(mv — a') ,S' 1 (my — a')> 



(27r)V(detE)l 



(27r)V(det E')s 



^ _ [<(mv — a),S 1 (mv-a)>-<(mv-o'),E l 1 (mv — a')>] 



/(my)l| mv |<s 



detE 



det £' 



(2.14) 



g— 5<(m\/- a),S 1 (mv— a)> 

(27r)V(detE)5 
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The last term is estimated by dropping the characteristic function and applying the first state- 
ment in Lemma 2.4. The last sup is estimated from above by 



1 - 



detE 
ditS 7 



det £ \ 2 
detX/y | m S ^|< s 



2 _ [<(mv— a),S 1 (m v —a)> — <(m v — a'),T,' 1 (m v — a') >] 



(2.15) 



where, using the simple estimate \e x — 1| < \x\e\ x \ =: g(\x\) the sup in the last expression can be 



estimated by g 



2 SU P|my|<S 



< (m^-a),S l {mv — a) > — < {my — a'), £' l (my — a') > 



For this last sup in the argument of g we use the upper estimates in terms of two-norms 
2 sup |< my,£ _1 (a - a') >| + |< a, £ _1 a > - < o^E'V >| 

\m v \<S 



+ sup 

|m v |<S 



< (my — a'), 



(my — a') > 



E" 1 -£'" 

< (25 + ||a|| 2 + |ja'|| 2 ) IIE" 1 ^ ||a - a'|| 2 + 2 (S 2 + ||a'||!) x HE" 1 - S ,_1 || 2 
Collecting our results gives Lemma 2.3. 

To apply the Lemma we just use the short notation 

a := o(A) := n v m™ 0A '^' hA = U v (1 - ^Aa)" 1 (m A (/i A ) + qd A , dA m dA ) 
a' := a'(A 2 ) := n r (1 - qA z .i)~ l (m^ 2 (/ia 2 ) , 0^d\ A J 



(2.16) 



(2.17) 



for the expectations of the |y (-dimensional Gaussians under consideration. We also denote by 
£ := Hy (1 — qA Z d) _1 Uy the infinite volume covariance matrix restricted to V, and corre- 
spondingly £' := E'(A) := U v (l - qA^^Uy. 

The volume difference of the covariances is fairly harmless: Given e > we can choose Aq 
sufficiently large, s.t. for all AdA , we have that 1 - ( de t xy (A) ) ^ — e an< ^ H S_1 ~~ S '(A) _1 1| 2 < 
e. Further, all matrix elements of E'(A) are bounded from above by the corresponding infinite 
volume expression E. In particular, we can use the upper bound TrE'(A) < TrE. (All of this 
can be explicitly seen from the random walk representation of the resolvent, see e.g. [K4]) 

Assuming these choices we get with Lemma 2.3 

J N [a; E] (dmy)f(my) - j M [a'; E'] (dm v )f(m v ) 



< 2 



\v\ 



e + (1 + e)g ((2S + ||a|| 2 + ||a'|| 2 ) W^h \\a - a% + 2e (S 2 + ||a'||l)) 



(2.18) 



(s-iMbr 



2 |1/| e~ " 2T ^" +2' y 'e 



(g-i^'ib) 2 



for ADA (e) where the E-terms appear now as fixed constants. 
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To estimate the ^-expectation of this bound we will decompose the space of the integer 
heights into a 'regular set' H := H(A 2 ,A) := H W (A 2 ,A) n H { - 2) {A 2 ,A) where 

H^(A 2 , A) := {h zd ; \\a(A)\\ 2 < B, \\a'(A 2 )\\ 2 < B} 
7i^(A 2 ,A) := {h z a- \\a(A) - a'(A 2 )\\ 2 < e 2 } 
and its complement. We get from this (with ||/||oo < 1) that 



(2.19) 



AT 



m 



rri aA ,u} A ,h A , 



; (1 - Q&a) (f)-M (1 - qA Z d)~ (m* A2 (/i A J , A c) ; (1 - qA Z d)~ (/) 



<^[W(A 2 ,A) c ] + 2 



1^1 



e + (1 + e) 5 ((25 + B + B) US" 1 1| 2 e 2 + 2e (S 2 + B 2 )) 



(S-B) 2 
g 2TrS 



(2.20) 



To estimate the exceptional set of integer heights we will prove below 
Lemma 2.5: 

(i) For all (arbitrarily small) 5 there exists a B < oo (sufficiently large) s.t. v A TL^ l \A 2l K) 
5 for all sufficiently large A 2 ,A. 



< 



(ii) For all (arbitrarily small) 5,e 2 , there exist choices of volumes A 2 cA (suff. large) s.t. 
v% [H 2 (A 2 ,A) C ] < 5 whenever AD A . 

But assuming this property, (2.20) can be made smaller than any given 5 for sufficiently 
large A in the following way: 

1) Choose B = S/2 large enough that a) the sum of the last two terms is smaller than 5/3 
and b) v A H^(A 2 , A)° < 5/6, according to Lemma 2.5(i). (So we must have that both A 2 , A 
are large enough.) 

2) Given these choices of S, B, choose e, e 2 small enough s.t. the middle term is smaller 
than 5/3. (Then the estimates hold true, after possibly enlarging A.) 

3) Finally there are then choices of A 2 cA s.t. v% H (2) (A 2 ,A) C < 5/6 for all AdA , 
according to Lemma 2.5(h). 

This finishes our discussion of the proof of Theorem 2; it remains however to give the Proof 
of Lemma 2.5: In fact, the Lemma holds under the following two weaker conditions of 



(a) Uniform integrability lim BToo sup A v% ]T eZ .i R 00 ; X , y \h y \ > B 







(b) Uniform summability lim RToo sup A J2 y ez d ;\y\>R R oo;x,y\h y \ > e 
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for each x G Z d . (These condition are implied from the hypothesis by Chebycheff, e.g. 



sup A v% 



J2y£Z d R<x>;x,y\h y \ >B< ^Y^yeZ* Roo;x,ySWP A l/%[\h y \].) 



(i) : Note that, due to the exponential decay of the resolvent we have for uniformly bounded 
boundary conditions that limA|oo Ily (1 — 9 a a) 1 qd A! d A rhdA = 0. It suffices to look at one 
matrix element, say a x (A) of the vector a(A). Using that Ra-,x, v < R<x>;x,y uniform in the 
volume, the form of m*(h), and the uniform boundedness of the random shift in the continuous 
spin Hamiltonian, it is immediate to see that what we need is implied by the above uniform 
integr ability condition. 

(ii) : The differences are estimated as follows. For x G V we have 



[(l-gA A ) 1 (m* A (h A ) + qdA,dAm dA ) ^ - (1 - qA Z d) 1 (m* A2 (h A2 ) , Z d\ A2 ) 

< E la- - a - i A ^)x!y\ K (M 



j/GA 2 



+ E (1 - <?A A )~* \m* y (h y )\ +qM (1 - ? a a)~* 



(2.21) 



< sup 

2/£A 2 



E a-^CK (Ml 



yeA 2 



+ E (1 - q^)x!y R (hy)\ + E (1 - Q^)xl 

y€A\A 2 y€dA 

We need to show that the ^-probability of the event that the r.h.s. is bigger than some £2 can be 
made small by choosing the volumes in a useful way. The last (deterministic) term converges to 
zero; so we assume that A is large enough s.t. it is smaller than ?2/3. Now, from (b) we know that, 
given any 5, we have for all sufficiently large A2 that Yl y eA\A 2 ~~ 1 A z d ) x \ \ m y (hy)\ — ^/3 
with (say) ^-probability bigger than 1 — 5/2. We fix such a A 2 . What we have just seen in the 
proof of (i) ensures that, for given S we can find a B such that the sum of y G A 2 on the r.h.s. 
of the last inequality is bounded by B, uniformly in A 2 , with (say) ^-probability bigger than 
1 — 5/2. Now it remains to choose A is large as we want to make the sup over y's in the fixed 
A2 as small as we want, and thus the first line on the r.h.s. smaller than £2/3 to finish the proof 
of Lemma 2.5.0 



Let us finally give the modifications needed to get the 

Proof of the Addition to Theorem 2:. Let us look again at first at the first term of the 
decomposition (2.8) where / is now a local observable that is only exponentially bounded. We 
introduce the same type of exceptional set H(A 2 , A). Then, after using the Lemma 2.4 for the 
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Gaussian expectation on the exceptional set, the analogue of (2.20) becomes 



rn 



A 

A 2 TrS 



;(l-?A£) '1 (/) -MUi- qA^y 1 (m\ 2 (/i Aa ) ,0 A =) ; (1 - gA^)" 1 ] (/) 



<2^e^^ 



e A||Q||2 l 1 



+ 2 \y\ 



\ R A 1 r 



M(A 2 ,A)' 

e + (1 + e)g ((2S + B + B) ^ + 2e (S 2 + £ 2 )) 



(S-B) z 
2Tr£ 



(2.22) 

To treat the ^-integral over the exceptional set, use the Schwartz inequality vf^ [e A " a " a 1h(a 2 ,a) c ] < 
(i/jl [e 2A H a ll 2 ] ) 2 Vj( [H(A 2 , A) c ] 2 . Now, given the assumption of the existence of exponential mo- 
ments, the first term on the r.h.s. can easily seen to be bounded by a constant independent of 
A. But after this, we are essentially in the same situation as after (2.20) and the way of choosing 
the parameters stays the same as before. 



We are still due the 



Proof of Lemma 2.1: By joint shift in height-direction we can assume that h x = 0. 

r — i _1 

We write the kernel in the form T x (h x = 0\m x ) = 1 + J2hez-,h^o fa( m x) with f h {m x ) := 
e (m-(h a; )-m;(^=o))m-i[m:(^) 2 -m:(h a; =o) 2 ]+ );iC (^)-^(^=o)_ Xo prove unicit y Q f the local max- 
imum of the kernel we note that m x t— ► Ylhez-h^o A( m aO (being a sum of strictly convex 
functions) is a strictly convex function; hence it has a unique local minimum which is the global 
minimum. 



To prove the bounds on the minimizer we look at the individual minimizers fh(h) of each 
of the pairs fh + f-h for h = 1, 2, — We will show that \m(h)\ < A, for all h. This implies that 
the minimizer of m x i— > Ylhez-h^o fh( m x) satisfies the same bound (since all terms in the sum 
are strictly decreasing [increasing] for m x > A [< —A].) Now, a computation gives 



m(h) 



log 



m* (2h + d x (h) - d x {-h)) I to [h + d x (-h) - d x (h = 0) 



h - d x (-h) + d x (h = 0) 



*\2 



+ 



(m*) 



[2h(d x (h) + d x {-h)) + d x (h) 2 - d x (-h) 2 ] - Vx (h) + r] x (-h) 



(2.23) 



Substituting the a-priori bounds of the random quantities \d x (h)\ < 5 d and \r] x (h)\ < S v we get 
from this 



\m(h)\ < 



1 



m*2h{\ - 5 d ) 



log 



+ { ^X[4hS d + S d ]+2S v 



< m* [46 d + <5 2 ] | 



4(1 -<5 d ) m*2(l-<5 d ) 




(2.24) 
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To see the last estimates, just look at the nominator Norm.(m x ) (see 1.4): It is simple to check 
that this sum converges and it is bounded from above as well as bounded from below away from 
0. From this the bounds (2.1) follow. <0 



III. A useful contour representation for discrete heights 

In this chapter we will treat the measures v A for the discrete height model that are given by 
(2.5) with (2.3). The effective finite volume Hamiltonian for the integer heights can be simply 
computed as a minimum of a quadratic expression in continuous variables: For any boundary 
condition m it reads 

inf iZ^A,u,A,ft A (mA) = -1 < m * A (h A ),R A ml(h A ) > A +\ ^2(m* x (h x )) 2 -J]^) 

Vd(Ac)(qm),R A m* A (h A ) > A 

1 q ^ 

- 2~ < fjd(Ac)(qrh),R A (f} 9 (A°)(qrh)) > A +- ^ fn 2 y with 

d(x,y) = l 

R A = (q- 1 -A A y 1 

(3.1) 

with ?7a(Ac)(m) := d A ^ A friQ A denoting the field created by the boundary condition. We also note 
that the continuous-spin minimizer is given by (2.6). 

We will deduce a contour representation for the v- measures. Let us give the commonly used 
notion of 'contour', adapted to this model: 

Definitions: A contour V in the volume A is a pair composed of a support TcA and a 

'height configuration' h A £ /Z A , such that the extended configuration (h A ,0 Z d\ A ) is constant on 
connected components of Z d \T_. A contour model representation for a probability measure 
v on the space Z K of integer height- configurations in A is a probability measure Q on the space 
of contours in A whose height-marginal reproduces v, i.e. v({h A }) = ^ r : Q({T}). The 

h A( r ) = h A 

connected components of a contour T are the contours 7, whose supports are the connected 
components 7 ofT and whose sign is determined by the requirement that it be the same as that 
of T on t~. 

The result of this chapter is 

Proposition 1: Suppose that q is sufficiently small, q(m*) 2 sufficiently large and 5d < \. 
Then there is a h A -independent constant K A (lo a ) s.t. we have the representation 



mt mA GIR A U A (m A ) 



K A (w A ) x e -<%).v(hA)> Y, Po(r>r) 



(3.2) 

h A( r ) = h A 
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for any h\ € Z, h . The quantities in the above representation are as follows: 

(i) Small fields: Sc is a random variable for each h 6 % and C<Z/Z d connected and we have 
used the notation 

< S(u),V(h A ) >■= Yl Yl Sc ^ where V h (h A ) := {x £ A,h x = h} ^ 
hezccAnv h 

The Sc(h) are functions of the random centerings dc(h) = (d x (h)) xe c for \C\ > 2. Up to a 
boundary term (see below), the single-site part (3S x {h) is the 'random depth' r\ x (K). 

We have the smallness properties, for all realizations of the disorder, 
\S c (h)\< Const sj e - constalcl , for \C\>3 with a = - log [1 + l/(2dg)] ~ - log - 

Zi — (j 

\S c (h)\ < Const (5 2 d + S v ) for \C\ < 2 

(3.4) 

with const , Const being of the order unity, depending only on the dimension. 

(ii) Contour- Activities: The activity po(T;u>r) is non-negative. It factorizes over the con- 
nected components ofT, i.e. we have pq(T;ujy) = Y\ 1 conn cp of r P° (t» > ) • 

For T not touching the boundary (i.e. dd\T_ = $) the value of po(T;uJr) is independent of A. 
We then have the 'infinite volume properties ' of invariance under joint lattice shifts of spins and 
random fields, as well as under joint shift in the height- direction. 

Peierls-type bounds: There exist positive constants (3,(3 s.t. we have the upper bounds 

0< Po (r>r) < e-^(r)-/3|£| (3 . 5) 
with the 'nearest-neighbor contour energy' 

E S (T) := ^2 \h x — h y \ for a contour Y = (v_, {Ji x ) x ehj 



{",s)cr 
\x — y\ = l 



In the above estimates the 'Peierls- constants' can be chosen like 



(3.6) 



d 



P = q(m ' )2 WT^f^FV ^C^x-njlogi,™-^! } (<a,0) (3.7) 
Probabilistic bounds for the small field: For \C\ = 1,2 we have IESc = and 

IP[\S c (h)\>t]<e~^ with a 2 = const [p\ + o\) (3-8) 
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Proof: To produce a sum of the type (3.2) we need to decompose the terms in the Hamiltonian 
in such a way as to exhibit a low-temperature part, a non-local field part and high-temperature 
parts that can be expanded. We will produce 'support of contours' from all these various 
sources. As we will see there we will have to introduce some type of support that occurs only 
for unbounded spin models with interactions having no finite range. 

Remember that we put zero boundary conditions. First of all it is now convenient to 
rewrite the integer-height Hamiltonian in the following form that makes explicit that it has 
purely ferromagnetic couplings: 



• n £ ^m 8A =0,UA,/lA 
m A £lR A A 



= ^2 Jk > x >v ~ h v~ dx(K) - dy(hy) ) - ^(^) + Kl (^ A ) with ( 3 - 9 ) 

{x,y}cA 

h x = h x l xe A, d x (h) = d x (h)l xe \ 



where Ki(u>a) is a constant that is independent of the height-configuration h&. The Ja's are 
positive and their nearest neighbor parts satisfy the lower bound for nearest neighbors and and 
upper bound for the decay of the form 



min J A;X;y >qm* 2 [4((l + 2dq) 2 -q 2 )] 1 

{x,y}<ZK,\x — y\ = l 



: ,, . (3-10) 



where \x — y\ is the 1-norm. Assuming this upper bound we have that 



Ja-,x, v < e a ' x y ' for \x — y\ > r(m*, q) where 



a = - log [1 + l/(2dg)] , r(m*,q):-- 



\og(m*\l + 2dq)/A) 
log(l + l/(2dg)) J 



+ 1 



(3.11) 



where square brackets means integer part in the definition of the integer range r = r(m*,q). 

Indeed, these couplings Jh;x,y can be conveniently read off from the following rewriting of 
the quadratic form appearing in the minimum of the continuous-spin Hamiltonian using the 
random walk representation of the resolvent: We decompose R\-x, y = SccA ^ ( x ~~ * V > ^) W1 th 
TZ(x — > y;C) = ^ 7 (g _1 + 2d) ( -' 7 ' +1 - ) where the sum is over all nearest neighbor paths 7 on 
the lattice from x to y that visit precisely the connected set C (see e.g. [K4]A.ll ff.) We have 
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then 

- < m A , R A m A > A +q ^ ' 



XI X ~^ y ' ,C ^2 [(m x - m y f - ml - m 2 y ^ +^ ^ ^ft(:r^y;C)m 

rr,yeA CcA " 16A CcZ d y£Z d 



x,y^A i£A 



X ft(x^y;C) 



C:Cn(Z< i \A)#0 



ml 



(3.12) 



Indeed, this gives immediately the closely related form 

inf H f 9K=0 ' WA > hA (m A ) 
m A eiR A 

= X Ja;x, v (h x -h y - [d x (h x ) - dy(hy)]f - ^2 VxiK) + ^2 K A . X (h x - d x (h x )f + K^ua) 

{x,y}cA x€A xGA 

(3.13) 

with 

jji*" 2 ^ m* 2 -\ 

JA-,x, y = ^- 22 TZ(x ^y;C), K ^ = ^~ 22 H{x^y-C) 

C:x,yeCcA C:Cn(z d \A)^9 

y.yGC 

To get the form (3.9) we can of course look at the boundary term as a coupling to a boundary 
condition h x = for x G <9A. Note that K A - X falls off exponentially as a function of the distance 
from x to <9A. Now, for every site x G A we pick a site y(x) G <9A that has minimal distance 
to x (with some arbitrary deterministic prescription to make this choice unique.) Then we 
extend the definition of J to all pairs in 2Z d x Z, d by Ja-,x, v '■= K A - X for y = y(x), J A -,x,y '■= for 
y G (A) c or x, y G (A) c . The lower bounds in (3.10) follow from the fact that 1Z (x — > x + e ; C) = 
[(<7 _1 + 2d) 2 — l] 1 for nearest neighbors x, x + e. The upper bound follows from the fact that 

/ 2d 

jtZ J v 7 

(see Appendix of [K4]). Now, for a given configuration fa a there are pair interaction terms 
JA;x,y {h x — hy — [d x (h x ) — d y (h y )]^ that are big (they will make up the essential contribu- 
tions to the 'low-temperature contours'), small (they will be expanded and make up high- 
temperature contours) and intermediate (they cannot be expanded and will be adjoined to the 
low-temperature contributions). The difficulty about these intermediate contributions is that we 
need to find conditions that ensure that their existence implies the existence of low-temperature 
contours nearby that will actually dominate them. 

Below we will introduce a set of pairs of 'big and intermediate interactions', £-^(h)cA x A. 
In particular the pairs that make up the low-temperature contributions will be contained in this 
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set of 'dangerous edges'. For fixed height configuration we decompose the exponential of (3.1): 

-J2 {x ,y } ^-(h) J^^yih^-hy-ld^h^-dyihy)]) 2 (3.16) 

E* 



e 



x e 



(E^cAr^ + E^Anv^ 4 I , 9 (iW-4C)) 2 + E, Wt 1'('') 

The rest of the proof is a careful treatment the three exponentials on the r.h.s. 



(i): First exponential in (3.16): Low temperature-contours (nearest neighbor parts, large 
fluctuation long range parts) 

Given h\, our aim is to define a support of a contour T^ T (h\), s.t. the first exponential 
can be written as a contour activity p LT (T LT (hA), h\) that satisfies a Peierls-type estimate in 
terms of the n.n. surface energy and the volume. We will have two contributions to the set of 
dangerous edges, £-^{h) ■= £^\h) U £^\h). For the first, the short range part we put 

s£\h) := {{x,y} G A x A;d(x,y) < r where h x ^ h y } (3.17) 

with the range r = r(q,m*) > 1 given above in (3.11). These interactions provide the 'main- 
part' of the low-temperature Peierls constant. Looking at these is in perfect analogy to what 
one would do in the case of an Ising model. The corresponding ingredient of the support of the 
LT-Peierls contours will be the connected components of the corresponding vertex set, i.e. 

(h) := {x G A; 3y G A s.t. d(x, y) < r where h x ^ h y } (3.18) 

Assuming that S < \ it is readily seen that we have an energetic suppression in terms of the 
n.n. surface energy: 

Ja;x, v \h x - h y - d x (h x ) - d y (h y ) J > m.n. ^ \h x - h y \ with 

{x,y}eS^\h) <x,y>eA (3.19) 

Tn.n. := (1 - 25 d ) 2 min J A , XjV 

{x,y}cA,\x-y\ = l 

Now we come to the second LT-part of the support of the contour, the large fluctuation 
long range part. It is unavoidable since the height variables can in principle have unbounded 
fluctuations: Indeed, if the difference between heights at far away sites is extremely large, it 
becomes impossible to treat the corresponding terms in the Hamiltonian by an high-temperature 
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expansion. It is however intuitively clear that such an event should be very unlikely since it 
implies a large energy cost due to the short range parts. To turn this into a contour representation 
we look at the following set of 'dangerous' bonds, 

£i 2) (h) := {{x, y}cA x A; d(x, y) > r; where \h x - h y \ > e% lx ~ yl } (3.20) 

whose interactions can not be treated by a high-temperature expansion, with a being the es- 
timate on the decay-rate of the J's, as given in (3.11). The corresponding part of the LT- 
support will have to contain the vertex sets of the connected components of the corresponding 
graph. Our aim is then to show the Peierls-type estimate in terms of the n.n. interface energy 
and the volume of the contour. The problem with the volume-estimate is that there is of 
coarse a gap between high-temperature and low-temperature expansions: if the interaction term 
Ja;x, v (h x — hy— d x (h x ) — dy(h y ) ^ is not small enough for a high-temperature expansion, 
the term itself need not be large enough to provide a low-temperature Peierls constant that 
is big enough. So, what might happen is that the large fluctuation long range parts just glue 
together connected components of the T^-parts without contributing much energy themselves. 
We will however show that in such a case the n.n. surface energy will be at least as large as 
the resulting volume of the total contour. Indeed, if the interaction term is not very small, each 
nearest neighbor path from x to y contributes a nearest neighbor interface energy of e%\ x ~ y \ . 
Patching together those paths along with safety cubes around them we will get the second part 
of the contour with a useful Peierls constant. The details are as follows. 

It is convenient to work in all of Z, d and define sets that will be the essential part of the low- 
temperature contours which are not necessarily subsets of A. The final supports of contours will 
then be obtained as intersections. Given our extended configuration h it is convenient to extend 
the above definition writing £ ( ^ d (h) := {{x,y}c2Z d xZ d ;d(x,y) > r; where \h x -h y \ > e a \ x ~ y \}. 
This is a finite set. To each pair {x, y} € £^\{h) we associate a cube Q ({x, y}) C/Z d among the 
cubes containing the points x and y with the smallest side-length. (If the line between x and y 
doesn't happen to be a diagonal, we choose some arbitrary deterministic tie-breaking procedure 
to make this choice unique.) Then we put T^ d {h) := \J^ X y y €£ (V^ Q y})- Finally we define 
the LT-support r LT (/i)cA of the contour by 

ir(h):=(T^\h)UTSl(h))nA (3.21) 



Then the desired lower bound of the n.n. surface energy in terms of the volume of the long-range 
parts is given by the following 

Lemma 3.1: Suppose that 7 is a connected component of the set T^ d (h). Then the nearest 
neighbor surface energy satisfies a Peierls estimate of the form 

e 2 \ /ea\ d 



£ \k -h y \>K |,| m th K := ^ f j > (_) (3.22) 



<x,y>Ej 
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Proof: From the set of cubes Q({x,y}) whose union makes up 7 we will consider in the 
following only the maximal ones w.r.t. inclusion. Let us denote them by Qi, i = 1,...,N. 
(That is we discard those that are contained in a strictly bigger one.) We prove the Lemma by 
induction over the number N of such cubes of a connected component. In the case of one cube, 
say Q(x,y) of side-length L, we have that 



^ \hi -hj\ > \h x - h y \ 

<i,j>eQ(x,y) 

> e fl*-vl > e § l > L d inf ( = 3 d K d L d 



(3.23) 



where the first inequality is the triangle inequality, the second the definition of the 'dangerous 
bonds'. This proves the single-cube case even with a constant 3 d K. 

Now, given a contour 7 = IJ^=i Qii we define a smaller contour r' in the following way. 
Pick one of the cubes Qi with the largest side-length, call this cube Qi a (with corresponding 
side- length L io ). Define 

T' := II Qi 

i€{l,...,iV};i^« \"-^J 

The resulting contour might be connected or not. Since we took away the biggest cube and its 
'contact partners' we have |r'| > | -y | — 3 d L d o . For the nearest neighbor surface energy we have 

\h x ~ h y \ > ^2 \hx ~ h y \ + \h x — h y \ ^ 25^) 

<x,y>&2 <x,y>€F' <x,y>€Qi 

Using the induction hypothesis on the connected components of T_' and the estimate on the 
single cube from above we get the desired 



\K -h y \ > K\T'\ + 3 d KL d o > K\j\ 

<x,y>€j 



(3.26) 







Now we are done with the treatment of the first exponential in (3.16): Indeed, the sum 



dx(h-x) — d y (hy) j l~ h _^~ h decomposes over connected compo- 



^2{x,y}e£-(h) ^ A ' X 'V y lx hy 

nents 7 of £ LT (/i). Denote the corresponding pairs by S-^.^h) (i.e. the pairs with vertices in 
7 U (7 n dh)). Denote by 7, the contour whose support is 7 and whose height configuration is 
h-£ on 77. We define the LT-activities by 



P LT (7*) := e 



(3.27) 



22 



The desired LT surface-energy /volume Peierls-type estimate is obvious: Denote T^ 2 - 1 := 7W n 
T^\h) the parts of the connected component that are due to the large-fluctuation- long-range 
part. Then we have 



^2 JA;x,y (hx -h y - d x (h x ) - dy(h y ) ) l' h ^ hy >r n , n E s (^i) 



-— EsM + S — (2r + l)' 



(3.28) 



>PE s { 1i ) + t 1 



7, 



with j3 = r " 3 "- and T\ = T " 3 "- x min{(2r + 1) d ,K}. Let us look at the large to*, small 
g-asymptotics with r„. n . ~ qm* 2 /4 sufficiently large, a ~ | log ^ r ~ 4 *° g "| . From this we get 

that /3 ~ qm* 2 /12, T\ ~ Const qm* 2 (^^^r^j ■ (Note that the contributions of the ET-term will 
actually obtain a better Peierls constant that hence won't be visible.) 

(ii): Second exponential in (3.16): High-temperature expansion 

Let us just write £ := £ LT (/i). We find it convenient to use a little rewriting of the exponent. 
Since we have exponential decay of the interaction, we can just 'fill the space between the 
endpoints' to define contours that obey Peierls estimates: To each pair {x, y} G [Ax A] \£-^(h) 
of sites we associate a 'one-dimensional' polymer g = g(x,y)cA that is the set of sites of one of 
the nearest-neighbor paths from x to y (with some prescription to make the choice of this path 
unique.) Then we have for the number of sites that \g\ = \x — y\\ + 1. We use this notation to 
denote the terms in the last sum by sets g and put 

Sg(hg) := l {XtyW£ _ {h) X J A . Xy y (h x -hy- d X (k X ) ~ dy (ky) ) H^ky ^ ^ 

We note that, due to the decay of the resolvent (with a), the uniform boundedness of \d x \ < Sd 
and the definition of the 'dangerous bonds' (with a/2!) we have < S g (h g ) < e - consta \a\. Note 
that there are only non-vanishing terms for g n £ / 0. (Indeed, in the case that both x and 
y are not in £, they must lie in different connected components of the complement of £.) The 
exponential can now be treated by the subtraction of bounds trick: We have 

*y conn.cp. of T 

(3.30) 

where n(£, g) counts the number of connected components of £ that are connected to g. The term 
under the first product defines a non- negative quantity r(-y) that is /i-independent and satisfies 
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1 > r{j) > e~'l' e (for suff. large a). The last exponential can be polymer-expanded 

and written as a sum Pr T1 (G, h G ) with a nonnegative activity < pr{G,hc) < 

e -const'a\G\ ( wn i cn [ s one f or empty G). 



(iii): Third exponential in (3.16): Non-local small fields 

The last exponential would not be present in the corresponding model without random- 
ness. It describes the random modulations of the 'vacuum-energy' caused by the ci-variables 
in the flat pieces outside the LT-contours (where height-fluctuations are occuring). To get the 
decomposition into Sc 's we use the decomposition of the resolvent and define for all C's 

l3S c {h c ) := -j- E K(x ^y;C)(d x (h x )-d y (h y )) ^ 3 ^ 

From the bound (3.15) we get 

(3Sc(h c )\ < ^) 2 ^?T^ q z- am - 2) (3-32) 

This is always fine for \C\ > 3; for \C\ = 2 we see that 5 d really needs to be small enough 
to get a useful bound. We put Sc '■= Sc for \C\ > 3 and Sc '■= Sc — IESc for \C\ = 2. 
[We remark that we could relax the assumption of smallness of 5d by the introduction of large 
d x (h) = d x (h)l\ dx ( h )\> Sl + d x (h)l\ dx ( h )\ <Sl =: d x (h) + d B x (h) . Then we would have to introduce a 
control field N x (h) that would contain a contribution of the type Const \d x (h)\ l| ( / a .(h)|>5 1 as well 
as a contribution from the large r/a,(/i)'s. Since this is simple but would obscure the structure of 
the contour-model we don't present the details here.] 

For a flat height configuration on C (i.e. if CcVh for an h G /Z) we just write Sc(h) := 
Sc(hc) the former defining the 'small field' appearing in the final contour-model representation. 
Then we have for the bulk term in the third exponential in (3.16) 

E E JA; X , y (d x (h)-d y (h)) 2 = (3j2 E S cW+(3 E S c (h c ) + K' A 

hez { x ,y}cAnV h heZCcAnVn cca v • > 

with some obvious ZiA-independent constant K' A . The first term is the desired small-field contri- 
bution. The second term is attached to the LT-contours and can be expanded. Using subtraction- 
of-bounds as before, its expansion gives 



cca Sc(h c ) 

Kq ^const 



II r 2 (j) 

7 conn.cp. of r LT (/l) 



E Pr™(h)(G,h G ) (3.34) 

G:Gnr LT (/ l )/0 
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with const a-decay. As far as for the bulk terms, the fields rj x (h) simply make up the local 
contributions to the small field. Finally, we discuss the corrections due to boundary effects. We 
write the interaction with the boundary in the form 

E E jA ' x <y ~ = - ^Vx(K) + K A where 

heZ Teel h X€A (3.35) 

fj x (h) := -K A . X (d x (hf -IE [d x (hf]) 

We note that |rfo(/i)| < qm* 2 ^p^f - The centered field fj x (h) will just give a small volume 
dependent modification of the local small field that we finally define by 

S x (h):= V l(h)+f) x (h) (3.36) 

We note that from this definition the probabilistic bounds for |C| = 1 are clear. Also, the 
probabilistic bounds for \C\ = 2 are obvious. 

Finally putting together our results from that (i)-(iii) we end up with the representation 
II (K7)r 2 (7)p LT (7,M) E P^ hA )(G,h G )p^ {hA) (G 2 ,h G2 ) 

7 conn.cp. of£ LT (ft A ) G: Gn£ LT ( h A ) ^0 

_ ' G 2 :G 2 nr LT (fc A )^0 

x e -<S,V(h A )> 

(3.37) 

Resumming over G, (SVs that have the same r LT (/i A )ur D (/i A )uGuG 2 =: r we get the desired 
form. This concludes the proof of the proposition. 



IV. Proof of Theorem 1 

The proof is a direct consequence of the representation of Theorem 2, the contour-representation 
for the discrete- height model of Proposition 1, and the results from the renormalization group 
analysis for a discrete- height contour model from [BoKl], [Kl]. It is crucial for this that the 
contour model constructed in Chapter III and given in Proposition 1 is 'renormalizable' with the 
procedure described in detail in [BoKl] in 'Chapter 4. The Gibbs State at Finite Temperature'. 
Indeed, it satisfies the inductive assumption of a contour model given in [BoKl] 4.1, p. 457, for 
the trivial choice of empty bad regions and vanishing control-field N. We must however have 
for this that the uniform bounds S v and Sd are sufficiently small; otherwise we would have had 
to introduce bad regions and large fields, which is however mainly a notational inconvenience. 
(There is further a completely trivial difference in that we have exponential decay for the small 
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fields only for |C| > 3; we could of course trivially cast the present contour representation into 
the one from [BoKl] by splitting the field Sc for C = {x, x + e} into new local small fields 
at sites x and x + e and producing a stochastic dependence up to distance 2.) The result of 
[BoKl] then gives that, for sufficiently large (3 (< a, (3), for sufficiently small of ff there exists 
a non-random subsequence of cubes A s.t. the measures v% (obtained from the zero boundary 
continuous Gibbs-measures) converge weakly to an infinite- volume Gibbs-measure for a.e. co 
[see [BoKl], p. 417, Theorem 1]. To conclude that convergence of the ^-measures takes place also 
on local observables / that are only polynomially bounded, \f(niv)\ < Const (1 + |my|) p , we 
would like to use the addition to Theorem 2 (2.7). Although its assumption on the convergence 
of the z/-measures on exponentially bounded observables is a very natural one that is believed 
to hold, it is unfortunately not a straightforward consequence of the RG-analysis. Along the 
lines of Chapter II, the reader will however have no difficulty to prove the analogous extension 
for polynomially bounded observables under the condition that sup A IEi/\\hy\ p < oo, for all 
exponents p. This assumption is in fact true; we even have (4.2), see below. 

Let us use the short notation M[h Z d) := TV (1 — qA Z d)~ 1 m* zd (h Z d) ; (1 — qA Z d)~ 1 
Then we have for the second moment 

JE^ «) = JEv w [M[h Z d] (m 2 X0 ) -N[h zd ] {m xo f}+IEv"N[h z a] (m XQ f 
= (1 - ?A^)-J + lEv" (1 - qA Z d)^ y m* y (h y )j 



(4.1) 



With the Schwartz inequality we bound the last expectation by (m*) 2 sup lEv 



(|/ij/l + fid) 



To estimate this last expectation, we utilize a corollary of the RG-analysis for the discrete height 
model of [BoKl] (whose complete proof can be found in [Kl], see Theorem 1.2 therein) saying 
that, for any q E /Z, we have 



JE^ [\h y \ q ] < K(q) (e + e - Const ^j (4.2) 



for any y, with g-dependent K(q) and some positive exponent k. This immediately proves 
(1.4).0 

To get more estimates on the continuous infinite volume measures in terms of the discrete 
one might now utilize ^[B] < ^[A ] + sup^ dGA M [h Z d] [B] for any good choice of an event 
A in integer height space. 
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